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Abstract. The existence of dwarf spheroidal galaxies with 
high internal velocity dispersions orbiting in the Milky Way 
raises questions about their dark matter content and lifetime. 
In this paper, we present an alternative solution to th e dark 
matter dominated satellites proposed by Ibata & Lewis ( 1998 ) 
for the Sagittarius dwarf galaxy. We performed simulations 
of two kinds of N-body satellites: the first models (f-models) 
could correspond to satellites with high dark matter content and 
they represent initially isolated models. The second models (s- 
models) have either low or negligible dark matter content and 
they are constructed in a tidal field. In spite of being on the 
same orbits, the s-models are able to produce a better agree- 
ment with some observational constraints concerning Sagittar- 
ius. From our simulations, we can also infer that Sagittarius is 
in the process of being disrupted. 

Key words: Galaxies: local group, interactions. Methods: nu- 
merical 



1. Introduction 

The current interest for dwarf spheroidal galaxies (hereafter 
DSphs) as satellites of the Milky Way is raised by two fun- 
damental questions concerning their evolution: the real content 
of dark matter (DM) in these objects (see, for example, Ma- 
teo 



1994,1997,1998; Piatek & Pry or 1995; Burkert 1997 and 



references therein) and their implication in the hierarchical for- 



mation of the galactic halo (e.g., Johnston et al. 1996). 

The Sagittarius (Sgr) dwarf galaxy is the closest known 
satellite galaxy to the center of the Milky Way, Rgc ~ 16 kpc 
(Ibata et al. 1995, 1997). Due to its proximity we can expect 
from its study an additional contribution to our understanding 
of DSphs in general. Since the announcement of its discovery 



by Ibata et al. (1994) the structure and evolution of Sgr have 



been extensively discussed and simulated by various authors: 



(1995), Whitelock et al. (1990, Mateo et al. (1996), Alard 



Ibata et al. (1995), Johnston et al. (1995), Velazquez & White 



( 1 996 ), Ibata et al. ( 1 997 ), Edelsohn & Elmegreen ( 1 997 ), Lay 



den & Saragedini ( |1997| ), Ibata & Lewis ( |1998[ ), Mateo et al. 
(1998). Important problems concerning this system are: 1) its 
possible lifetime before dissolution, 2) the possible presence of 



Table 1. Observational parameters of Sagittarius dwarf galaxy 



(Ibata etal. 1997) 



Parameter 


rhb 


0.55 kpc 


o-o 


11.4 km/s 


fJ-oV 


25.4 mag/arcsec^ 


Lt 


> 10^ L0 


{M/L)o 


5OM0/L0 


Mt 


> 10^ Mq 




(5.6°, -14°) 


{u,v,wj^ 


(232,0:, 194)± 60 km/s 




25 kpc 


Rgc^ 


16 kpc 




171 kms-i 


{dv/dbj^ 


-3 km s~^/degree 



" Galactic coordinates 

Galactic velocities 

Heliocentric distance 
'' Galactocentric distance 

Radial velocity 
^ Gradient of the radial velocity 



DM in it. A critical point in the studies mentioned above is the 
question whether Sgr is in virial equilibrium or not. 

Usually, the total inferred mass of the DSph is calculated 
by assuming it is in virial equilibrium!. ^'^ ^^^^ case, the central 
mass-to-light ratio depends on the velocity dispersion through 
the equation (Richstone & Tremaine |l986| ) 



9ai 



2TTG^j.orhb 



(1) 



where 77 is near unity for a wide variety of models, (Tq the cen- 
tral velocity dispersion, fig the central surface brightness and 
Thb the half-brightness radius. The analysis of the validity of 
Eq. (|l]) for evolving DSphs has been studied in detail in the 
present context by Kroupa (1997, hereafter K97). 

In Table |l], we have summarized the parameters of Sgr 
DSph measured by Ibata et al. (1997). For the values of (Tq, 



In this paper, the inferred mass is always obtained by assuming a 
virial equilibrium for the satellite. 
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Thb and fioV given in Table |l], those authors obtain a high cen- 
tral mass-to-luminosity ratio {M/L)o — 50 Mq/L© by using 
the Eq. ([l|). Therefore, the total inferred mass assuming virial 
equilibrium, Mt, is DM dominated. The values of total lumi- 
nosity, Lt, and the total inferred mass, Mt, are also listed in 
Table |l]. However, the values of {M/L)o and Lt suggested by 
Mateo et al. (1998) are different. These authors have discov- 
ered a tidal extension of Sgr, which implies Lt < 5.8 x 10^ 
L0 and, as a consequence, the {M/L)o ratio could decrease to 
10 Mq/Lq. Adopting the structural and orbital chara cteristics 
usually assumed (given in Table |]), Ibata et al. ( 1997 ) recently 



suggested that such a satellite galaxy is expected to be tidally 
disrupted and destroyed after several pericentric passages, un- 
less a significant quantity of DM is present inside it, as inferred 
assuming virial equilibrium. The observations of the Sgr glob- 
ular clusters give an age spread of its constituents > 4 Gyrs 
and the youngest globular cluster (Terzan 7) is 9-13 Gyrs old 
(Montegriffo et al. 1998). The fast disruption obtained by the 



numerical methods raise the question about the age and the dy- 
namical history of Sgr However, it must be emphasized that 
the initial time t = of the simulations does not necessarily 
coincide with the time of formation of the oldest constituents 
of the satellite. 

The partial formation of the galactic halo by hierarchical 
processes (accretion of small galaxies) has recently received an 
increase of interest, stimulated for instance by the investigation 
by Lynden-B ell & L ynden-Bell (1995) on the reality of streams 
(Lynden-Bell |l982| ) in the close environment of the Milky Way, 
among them the well known Magellanic Stream. DSphs seem 
to belong to one or another of these streams and their evolution 
clearly depends on their environment. This scenario of satellite 
formation involves a low DM content for the DSphs (Barnes & 
Hernquist 1992; Kroupa 1998). However, the only simulations 
of no-dark matter satellite galaxies able to survive in a tidal 
field are those of K97 and Klessen & Kroupa (1998, hereafter 



KK98). These authors have studied a region in the parameter 
space (Msatj'Tsat), where Msat and rgat are the mass and a 
typical radius of the satellite, and they have obtained a residual 
satellite from a more massive one. Nevertheless, their remnant 
satellite galaxies are fainter than a typical DSph, unless the true 

(M/L)real < 3. 

These circumstances make questionable the maintenance of 
dynamical equilibrium and consequently the virial estimation 
of mass for these systems. 

In this paper we present self-consistent N-body simulations 
of Sgr and a more accurate and plausible scenario of its evolu- 
tion, which suggests that Sgr is likely to be able to survive for a 
long time (6 — 10 Gyr) by orbiting in the Galaxy without being 
dominated by DM. 

In Sect. 2, we present the model of the Milky Way used 
in our simulations. In Sect. 3, we describe the different mod- 
els of the satellite galaxy and corresponding scenarios of inter- 
action with the Galaxy. In Sect. 4 some numerical considera- 
tions are given. In Sect. 5, we present the results of our sim- 
ulations for the different chosen scenarios. These results are 
discussed in Sect. 6 in connection with recent conjectures by 




R [kpcj 



Fig. 1. Rotation curve of the initial Milky Way model (full line), 
with the contributions of each component 



Kroupa ( 1998b , hereafter K98) on the parameter space possi- 
ble for progenitors of surviving DSphs without DM. Finally, a 
summary is given in Sect. 7. 

2. The model of the Milky Way 

The initial model of the primary galaxy, representing the 
Milky Way, is based on the axisymmetric initial conditions 
adopted by Fux (1997) in his N-body modeling of the Galaxy. 
The mass distribution is divided in three components: (i) an 
oblate stellar nucleus-spheroid with a spatial density p oc r^^-^ 
in the central region, in agreement with the near-infrared ob- 



servations of the inner bulge (Becklin & Neugebauer 1968 
Matsumoto et al. 1982), and oc 7-^3.3 outside the bulge, as the 



number density counts of stellar halo objects (e.g. Preston et 
al. 1991 for RR Lyrae; Zinn 1985 for globular clusters), (ii) a 



double exponential stellar disc, with radial and vertical scale 
lengths hji = 2.5 kpc and = 250 pc and (iii) an oblate dark 
halo with an exponential profile ensuring a roughly flat rota- 
tion curve out to i? = 40 kpc. Both oblate components have 
a flattening c/a — 0.5. The total mass of the luminous com- 
ponent is Ml = 8.25 x 10^° M©. The model is identical to 
Fux's model m04t0000, except that the dark halo is more ex- 
tended, with a scale length 6=13 kpc and a total (untruncated) 
mass Mdh — 2.4 x 10^^ M©, and that the truncation radius 
is moved further away at Rc = 70 kpc. The chosen truncation 
radius ensures a non-vanishing density everywhere along the 
orbit of the satellite galaxy. The resulting rotation curve of the 
model is shown in Fig. |l]. 

The initial kinematics is obtained by solving the hydrody- 
namical Jeans equations, assuming an isotropic velocity disper- 
sion for the disc component and a anisotropic but centrally ori- 
ented velocity ellipsoid for the other components. The velocity 
distribution of the nucleus-spheroid and dark halo is generated 
from a 3D generalization of the Beta distribution which limits 
the number of escaped particles and thus greatly improves the 
equilibrium of the isolated dynamical model. More details are 
given in Fux (1997). 
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3. Satellite models 

The effects that a DSph suffers when it is accreted on a primary 
galaxy strongly depend on its structure and its orbit. 

We have considered two kinds of satellite models. They 
correspond to two scenarios of interaction between the Galaxy 
and the satellite. 

In the first models (f-models), we represent the satellite as 
a standard King's sphere (King 1966), which matches the ob- 
servational constraints for the total virial mass, Mt, the core 
radius. To, and the central velocity dispersion, do, of the Sgr 
DSph. In this case, the effects on the DSph are maximal be- 
cause the satellite, originally in an isolated situation, suddenly 
undergoes strong tidal perturbations. 

Our alternative scenario (s-models) assumes that either the 
DSph is formed in the tidal tail of another major accretion event 
and therefore it is built in equilibrium with the environment, or 
the DSph falls slowly from a quasi-isolated situation to a tidal 
field region and it has enough time to readjust itself to the envi- 
ronmental forces. In both situations, we begin our simulations 
when the satellite has already reached the equilibrium with the 
dense environment and, therefore, the life time of the DSph 
orbiting in the galaxy is expected to be longer than in our first 
scenario. In the present case, the satellite is modelled by a mod- 
ified King's sphere (Gomez-Flechoso & Dommguez-Tenreiro 



1998), which takes into account the tidal potential produced by 
the primary galaxy. 

We have checked several orbits in order to compare the ob- 
servational features of the Sagittarius dwarf and the numerical 
results in both fast and slow accreting scenarios. The orbits we 
have chosen reproduce the present position, {I, b), and galac- 
tocentric velocity, {U, V, W), of Sagittarius DSph in the obser- 
vational range (Table |lj), and, therefore, the simulated orbits 
are polar, as suggested by the void component of the proper 
motion parallel to the Galactic Plane. The orbits are either low 
eccentricity orbits in the central part of the Galaxy or high ec- 
centricity orbits. 

The central mass-to luminosity ratio of the models has been 
calculated by using the Eq. ([l|). The central velocity dispersion 
(To involved in this equation has been measured along the line- 
of-sight. We have removed the particles with the largest radial 
velocities (relative to the radial velocity of the center of mass of 
the satellite) to prevent contamination by outlying particles. We 
have only considered those particles with projected distance to 
the center of the satellite smaller than 0.5 kpc. 



Table 2. Physical parameters of the isolated King's model of 
Sgr (f-models): core radius, Tq, central velocity dispersion, ctq, 
total mass, Mt, dimensionless potential. Wo and tidal radius, 

n 



To 




Mt 


Wo 


rt 


(kpc) 


(km/s) 


(10^ Mo) 




(kpc) 


0.527 


15.0 


12.0 


3.26 


2.736 



Table 3. Parameters of the orbits of the f-models 



Model 






Period 




(kpc) 


(kpc) 


(Gyr) 


f-A 


12 


18 


0.23 


f-Bl 


8 


38 


0.45 


f-B2 


10 


70 


0.95 



ulation. Once in the inner orbit, the satelUte will evolve quickly, 
because tidal forces are strong in these regions. This scenario 
is unrealistic, because, in reality, a satellite which reaches the 
denser parts of the galaxy has suffered the influence of the 
galaxy potential for a long period of time. However, we have 
run simulations in such a case because they correspond to the 
common initial conditions assumed in the literature (Johnston 
et al. 



1995 



Oh et al. 1995 Velazquez & White 1995, Johnston 



et al. [19961 Edelsohn & Elmegreen [1997) K97, KK98, Ibata & 
Lewis 1998) and in order to compare the results with those of 
the more realistic s-model satellites. 

The parameters of the King's model we have selected for 
the satellite model in a fast accreting scenario are given in 
Table || and they have been chosen to reproduce the present 
characteristics of Sagitta rius, a s was done in other simula- 
tions (Velazquez & White 1995 ). The total mass of the satellite 
model corresponds to the virial mass inferred for Sgr There- 
fore, if the observational constraints on the luminous mass are 
considered, the satellite model could represent a DM domi- 
nated DSph. For these models, we have assumed {M/ i)]-eal = 
10 M0/L0 for the calculations of surface brightness, which is 
the lower limit of the mass-to-luminosity ratio for DM dom- 
inated satellites. A higher {M/L)^.^^^ would represent fainter 
satellite galaxies (for the same DM content) and an apparently 
faster dissolution process. 

The initial apocenter and pericenter and the period of the 
low eccentricity orbit (f-A) and the high eccentricity orbits (f- 
Bl and f-B2) are fisted in Table |. 



5.7. f-models 



3.2. s-models 



The distribution function of an isolated galaxy fulfills the col- 
lisionless Boltzmann equation. We can represent an isolated 
DSph as a solution of this equation. King's spheres are an ex- 
ample of that. 

If an initially isolated DSph reaches the inner regions of a 
galaxy within a short timescale, it has no time to modify its in- 
ternal structure. In this case the satellite maintains its isolated 
King's sphere distribution function at the beginning of the sim- 



In the other possible scenario we assume that either the satellite 
has been formed inside the tidal tail of a major merger (e.g. 
numerical simulations by Barnes & Hernquis t 1992 , and the 
observational counterpart by Due & Mirabel 1997 ) or it has 
been slowly accreted. 

In the first case, if the dwarf galaxy is formed in equilibrium 
with the tidal force of the environment, it does not contain a 



significant amount of DM (Barnes & Hernquist [1992[ ). 
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Table 4. Galaxy parameters of the s-model satellites 



Model 


ro 


(Jo 


Mt Wo 


rt 




(kpc) 


(km/s) 


(10' M©) 


(kpc) 


s-A 


0.06 


11.0 


0.93 3.78 


0.55 


s-Bl 


0.1 


11.0 


1.66 3.95 


0.99 


s-B2a 


0.1 


15.0 


5.27 5.42 


1.93 


s-B2b 


0.1 


11.0 


2.33 4.87 


1.48 


s-B2c 


0.3 


15.0 


6.04 2.98 


2.02 


Table 5. Orbital parameters of the s-model satellites 


Model 




fmax 


Tave Period 






(kpc) 


(kpc) 


(kpc) (Gyr) 




s-A 


15 


20 


20 0.25 




s-Bl 


8 


40 


35 0.45 




S-B2 


10 


74 


50 1.25 





In the second case, according to cosmological models of 
hierarchical structure formation, satellite systems are produced 
around massive gala xies. T hese satellites could contain DM ha- 
los (e.g. Cole et al. 1994 ). Tidal forces could be negligible in 
the outer regions of the galaxy where the satellite is formed. 
However, if we assume that the satellite does not go through 
the denser and central regions of the main galaxy in the first 
perigalacticon, we could expect that an initially massive satel- 
lite falls slowly on the center of the galaxy by dynamical fric- 
tion, loses part of its mass and reaches a central orbit. In this 
way, the satellite has time to modify its internal structure and 
to reach the equilibrium with the environment. 

DSphs described in the two last scenarios have to be in 
equilibrium with the tidal forces of the environment. In the 



paper by Gomez-Flechoso & Domfnguez-Tenreiro (1998), the 
structural parameters (total mass, Mt, velocity dispersion, ctq, 
core radius, Tq, etc.) of a satellite in the tidal field of the primary 
galaxy have been estimated. Those authors have proved that, in 
general, a galaxy in a tidal field can be described by a two pa- 
rameter distribution function. They have solved altogether the 
Poisson equation and the collisionless Boltzmann equation for 
a galaxy, taking into account the potential of the galaxy and 
the tidal potential of the environment. Only spherical terms of 
the tidal field were considered in this theory. However, the ob- 
tained equilibrium solutions are better representations of the 
system than isolated models. 

This result suggests that in our problem we could repre- 
sent a DSph in equilibrium with the tidal field of the primary 
galaxy as a modified King's sphere (see Gomez-Flechoso & 



Dommguez-Tenreiro 1998) with two free parameters. We have 
chosen the central velocity dispersion, (Tq, and the core radius, 
To, as input parameters, because both can be determined from 
observations. In our simulations, the values for these two pa- 
rameters are Tq ^ 0.06 — 0.3 kpc and CTo ^ 11 — 15 km/s, 
which reproduce the characteristics of the Milky Way satellites. 
Thus, we will try to explain the Sagittarius satellite as a typical 



DSph which has evolved orbiting for a long time in the Galaxy 
potential. The other parameters of the model (total mass, tidal 
radius and dimensionless central potential) are automatically 
determined by the tidal potential at each position of the orbit. 

The satellite parameters of the models which have been per- 
formed are listed in Table The second column is the core ra- 
dius, To, and the third column is the central velocity dispersion, 
ctq, which are input parameters of the modified King's spheres. 
The total mass, Mt, the dimensionless central potential. Wo, 
and the tidal radius, rt, (columns 4, 5 and 6) are output param- 
eters obtained by solving the collisionless Boltzmann equation 
with the tidal potential at the averaged distance of the orbit to 
the Galaxy center (parameter rave in Table ||). 

As it can be seen in Table ^ the dimensionless central po- 
tential and the mass of the satellite decrease for inner positions 
of the equilibrium satellite for models with the same Tq and (Tq. 
Furthermore, the mass of the models is smaller than the mass 
inferred from observations using kinematic arguments (Ibata 
et al. 1997, Mateo 1994) and it is in agreement with the ob- 
served luminous mass, assuming {M/ L)^^^^ ^2 — 5 M© /Lq. 
Therefore, we have assumed (-^/i)real ^ ^ ^o/Lq for all 
the s-models, that is a typical value for the stellar population 
of a DSph. The low value of the mass-to-luminosity ratio is in 
agreement with satellites either formed in tidal tails of major 
accretion events or tidally modified by orbiting for a long time 
in a tidal potential. 

The apocenter, rmax, the pericenter, rmin, and the period 
of the orbits are listed in Table ^. The s-A orbit is an orbit of 
low eccentricity in the inner region of the Galaxy and s-B 1 and 
S-B2 orbits have higher eccentricity. 

4. Numerical details 

The models have been evolved us ing th e treecode algorithm 
kindly provided by Barnes & Hut ( 1986 ) with a tolerance pa- 
rameter 9 = 0.7 and a time-step At = 1 Myr. The number of 
particles of the luminous and dark halo components of the pri- 
mary galaxy are 15671 and 29648, respectively, and the mass 
of the dark matter particles is 3 times larger that for the lumi- 
nous particles. All the satellite models have 4000 equal-mass 
particles, except s-B2a and s-B2c which have 8000 particles. 

We use a softening length varying proportionally to the cu- 
bic root of the particle mass of the component, in order to 
avoid well-known usual numerical effects in the simulations 
(e.g. Merritt 1996| ; Theis & Spurzem 1999). For the luminous 
particles of the Galaxy it is eL = 0.23 kpc and for the dark 
matter halo edh = 0.33 kpc. The f-models have £5 = 0.06 
kpc, but this value is changed to es — 0.05 kpc for the s-A, 
s-Bl and s-B2b satellites and es = 0.04 kpc for the s-B2a and 
s-B2c satellites. 

5. Results 

5.7. The main galaxy 

The main galaxy develops a bar-like structure, described by 
Fux (1997) for an isolated model of our Galaxy. 
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The global effects of the satellite on the primary galaxy are 
weak, since the mass ratio of both objects is huge. Besides, 
the poor resolution of the Galaxy model prevents a detailed 
description of the local effects of the Sagittarius accretion on 
the Milky Way. Therefore, we only deal with dynamical effects 
felt by the satellite galaxy. 

5.2. f-models 

For our f-models, we begin the simulations when the satellite 
has already reached the inner regions of the primary galaxy. A 
DSph galaxy on a low eccentricity orbit in these inner dense 
regions of the Galaxy undergoes stronger disruptions than on 
more eccentric orbits, because the tidal forces are stronger at 
all positions on the trajectory. 

5.2.1. f-A satellite 

In Fig. ^a (b), we plot the angular distribution of the f-A satel- 
lite along (perpendicular to) the orbit (which has eccentricity 
e = 0.2) as seen from 24 kpc away (it corresponds to the So- 
lar neighbourhood viewpoint), for four snapshots. In this case, 
the lifetime of the dwarf galaxy, before significant disruption, 
is nearly 0.4 Gyrs (Fig. ||a). 

The material from the satellite galaxy is tidally stripped, 
forming a long stream and then, after 1.2 Gyrs, an almost close 
great circle in an Aitoff projection, which survives for a long 
time (Fig. ^). Moreover, as we can see in Fig. the mean 
width perpendicular to the orbit on the sky is 7° (it was only 
1.9° at the initial time) and the projected surface brightness 
(Fig. ^) is 5 mag fainter at the end of the simulation (after 2. 1 
Gyrs). 

5.2.2. f-Bl andf-B2 satelHtes 

Forbothf-Bl and f-B2 orbits (e = 0.64 and 0.75, respectively), 
disruption mainly occurs at perigalacticon, because tidal forces 
are more efficient at small galactocentric radii. The lifetime of a 
satellite depends strongly on the orbit shape. The f-Bl satellite 
is destroyed after 0.5 Gyr (Fig. ^) whereas the f-B2 satellite 
survives for 1 Gyr (Fig. ^p) due to the longer period of its orbit, 
at this time the f-B2 satellite undergoes a close interaction (r ~ 
2 kpc) with the center of the primary galaxy and it is tidally 
destroyed. 

The final destruction of dwarf galaxies is more efficient 
in our more eccentric orbits and they present fainter surface 
brightness at the end of the simulation than the satellite on a 
low eccentric orbit (Fig. |]). The satellite particles of f-B mod- 
els are spread over all directions and no predominant streams 
are formed, as it can be seen in Fig. ^ for the f-Bl model. In 
Figs. |d and |f, we have represented the width (perpendicular 
to the orbit) of the stream. At the end of the simulations, the 
f-B2 models do not have any predominant peak in the mass 
distribution perpendicularly to the initial orbit. 

As a general result, due to the tidal field on the satellite we 
observe: i) a modification of the internal structure of the satel- 



I f-A model 



Gyr '-- model 

---- 0,315 Gyr:: 

- - 0.4^ Gyr -: 

■ -■ 3.1 Gyr " 



Gyr 

---- 0,315 Gvi 
- - 0.4^ Gyr 
■ -■ 2.1 Gyr 



c) f Bl modelf 



Gvr 
0.42 Gvr 
0,525 Gyr:: 
2.1 Gyr 



d) f HI mode 



Gyr 

---- 0.42 Gyr 
- - 0.525 Gvr 
--- 2.1 Gyr 



) f-B2 modef^ ^y^, - f) f-B2 model T 



1945 Gyr' 
- - 1.05 Gyr 
.- 2.1 Gyr " 



-- - 0.945 Gyi 
- - 1.05 Gyr 
2.1 Gyr 



-15 10 5 D 10 15 15 10 5 5 10 15 

(degrees) 

Fig. 2. f-models: mass distribution parallel (a, c and e) and per- 
pendicular (b, d and f) to the orbit as function of the angle (j) 
to the satellite center, as seen from 24 kpc away, which is the 
present distance between the Sun and the Sagittarius DSph 

lite for both cases of orbits (eccentric and quasi-circular), with 
increase of the projected velocity dispersion, and ii) a contin- 
uous loss of satellite mass and luminosity. Anyone tempted to 
infer a (M/ L)o ratio from such experiments by assuming virial 
conditions would find values 100 — 200 times higher than the 
true ones at the disruption time (Fig. ||), confirming the results 
ofK97. 

5.3. s-models 

In this subsection, we analyze the interaction effects on a DSph 
in equilibrium with the galaxy potential. Either it could have 
fallen down slowly from an intermediate region in the denser 
parts of the primary galaxy, loosing part of its mass and becom- 
ing a low DM satellite, or it could have been formed in a tidal 
tail of a major accretion event. 

5.3.1. s-A model 

DSphs which are theoretical equilibrium solutions to Sgr at 
a circular inner orbit are small and low mass galaxies. The 
tidal field is almost constant along this quasi-circular orbit and, 
therefore, interaction effects on the satellites in equilibrium are 
not important. For the s-A orbit, there is a continuous loss of 
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f-A model {2.1 Gyr.) 
b=90° 




t-Dl model (2,1 Cyr,) 
b=90° 














90° 1^ '^180° ■ f 






870° ^ J 














b = - 


90° 



t (Gyr) 



Fig. 4. Aitoff projection of the f-A and f-Bl models at the end 
Fig. 3. Radii which enclose 70%, 60%, 50%, 40%, 30%, 20% °^ simulation (2. 1 Gyr) 
and 10% of the initial mass of the satellite, for f-models 



mass, but the satellite has still ^ 20% of the initial mass (Fig. 
^) after 2 Gyrs and it is still detectable (Fig. ^a). The cen- 
tral surface brightness has only changed by 2.5 mag (Fig. ^a), 
evolving from 21.3 mag arcsec"^ to 24.8 mag arcsec"^. The 
line-of-sight velocity dispersion grows in the outer parts of the 
satellite in agreement with findings by K97, remaining almost 
constant in average in the inner parts (< 0.2 — 0.3 kpc) (Fig. 
^). We have calculated the mass inferred from the line-of-sight 
velocity dispersion, measured inside a radius of 0.5 kpc from 
the satellite center, using Eq. ([l|). Since the central velocity dis- 
persion and the projected density (Fig. ^) evolve slowly in the 
central region of the satellite at the beginning of the simulation, 
the inferred mass-to-luminosity ratio remains almost constant, 
contrary to the f-models. However, as we have explained above, 
the real mass of the satellite decreases and it is only ~ 20% of 
the initial value after 2 Gyrs. The structural evolution is slow 
(Fig. ^) and it leads to a lower central surface brightness at 
the end of the simulation and a {M / L)o ratio calculated from 
the velocity dispersion (inside 0.5 kpc) which is 5 times higher 
than the real one (Fig. ^). This effect could be dramatically 
increased if the velocity dispersion is measured inside a larger 
radius, because of the velocity dispersion increase in the outer 
region of the satellite, leading to a calculated {M/L)o ratio up 
to 10 times larger than the real one. Very faint tidal streams are 
formed, which are spread along ^ 75 — 100% of the orbit with 
a spatial width ^ 6° (Fig. ||). 




Gyr. -- 
0,525 Gyr 
3,1 Gyr, -- 



Gyr, 
1,05 Gyr 
2 1 Gyr, 




-0,5 0,5 -1 -0,5 0,5 -1 -0,5 0, 

log(r) (kpc) 

Fig. 5. Surface brightness /iy of the f-models, assuming 
(W^)realy ^lOMo/Lo 



5.3.2. s-B models 

Perigalacticon and apogalacticon of high eccentric s-B 1 and 
S-B2 orbits decrease with the time. The most dramatic exam- 
ple is the s-Bl orbit. The DSph passes 1 kpc from the galaxy 
center at its third passage at perigalacticon (Fig. ^), because 
of the orbital energy loss. During this passage the satellite is 
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Fig. 6. Evolution of the log of the ratio between the mass-to- 
luminosity relation, {M/L)o, inferred from the velocity dis- 
persion values under dynamical equilibrium conditions, and the 
real mass-to-luminosity relation, (A//i)j.g^j, for the f-models 



strongly tidally stripped and it loses most of its mass (Fig. ^). 
The surface brightness of the satellite decreases with time (Fig. 
^) and the velocity dispersion and the mass-to-luminosity ra- 
tio increase in the outer parts (Figs. ^ and |^). The satellite is 
finally completely destroyed after nearly 1 Gyr 

S-B2 is the most external and eccentric orbit which we have 
chosen. As the other models, it fulfills at certain times the po- 
sition and velocity constraints for t he Sag ittarius DSph orbit in 
the observational range (Ibata et al. 1997). 

In order to investigate this case in more details, we have 
selected three different satellite models for orbiting on this tra- 
jectory (s-B2a, s-B2b and s-B2c). The concentration (defined 
as c = \og{rt/ro)) varies from 1.28 (s-B2a) to 0.83 (s-B2c). 
The concentration determines the fate of the satellite, as well 
as the evolution of the {M/L)o ratio. The more concentrated 
the satellite is, the less variation of the inferred {M/L)o ratio it 
suffers (Fig. 

The {M/L)o ratio in Fig. |l^ has been calculated for the 
central region of the satellite (< 0.5 kpc), where the tidal ef- 
fects on the mass and the velocity dispersion are weaker How- 
ever, even in the most concentrated satellite (s-B2a), we obtain 
{M/L)o ratios which are 10(A//L)j.g^j, these values rise up 
to several hundreds for the least concentrated model (s-B2c). 

It is interesting to notice the strong variations of the inferred 
(M/L)o ratio with time. This behaviour is caused by the evo- 
lution of the surface mass density when the satellite suffers a 
close encounter with the core of the primary galaxy. It leads 
to a successive rearrangement of the internal structure of the 
satellite. Another important parameter in the {M/L)o calcu- 
lations is the angle between the observer-satellite line and the 
main axis of the satellite. The strong anisotropics of the satellite 
(tails along the orbit, anisotropic velocity dispersion, etc) could 
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Fig. 7. s-A model: a surface brightness assuming 
(Af/L)j.g^j y = 2 M0/L0, b velocity dispersion of the satel- 
lite as a function of the radius as same timesteps as in a), c 
evolution of the radii which enclose from 70 to 10 % of mass 
of the satellite from the top to the bottom, and d evolution of 
the mass-to-luminosity ratio measured from the velocity dis- 
persion. The observational values of /i, a are given in Table 



produce various {M/L)o values (as already suggested by K97 
and KK98). 

The line-of-sight velocity dispersion evolves as shown in 
Fig. [ll]. The projected velocity dispersion decreases in the cen- 
tral part of the satellite (r < 0.5 kpc) and it increases for 
r > 0.5 kpc. Thus, the (M/ L)o ratio calculated from Eq. (|l|) 
varies, depending on the limit radius used for the measurement 
of the central velocity dispersion (this limit radius is related to 
the observational resolution). 

s-B2a model: 

The s-B2a model undergoes a long term evolution. It is 
the most concentrated model that we have chosen and it sur- 
vives for at least 10 Gyrs. A mild loss of mass is produced, 
mainly at perigalacticon: the satellite looses the outer layers of 
mass, but a core enclosing 30% of the initial mass subsists (Fig. 
p^). The limit radius is smaller than the initial one, but the 
half-brightness radius increases (the final satellite is less con- 
centrated). The central surface brightness decreases (Fig. [l3|a). 
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Fig. 8. a Aitoff projection of the s-A model at the end of the 
simulation (2. 1 Gyr), and mass distribution b parallel and c per- 
pendicular to the orbit as function of the angle (f) to the satellite 
center 



changing from ~ 21.0 mag arcsec^ to ^ 22.0 mag arcsec^ 
(assuming (M/L)j.g^j = 2 Mq/L©). The final s-B2a satel- 
lite is smaller {rhb — 0.13 kpc) than the observed Sgr DSph 
{fhb = 0.55 kpc) and the {M / L)o ratio is not large enough to 
reproduce the inferred {M/L)o — 50 Mq/L© by Ibata et al. 
(|1997|). 



s-B2b model: 

The mass loss of the s-B2b satellite (Fig. p^) is stronger 
than that of the s-B2a satellite. The particles of the satellite 
outer region are stripped, mainly at perigalacticon. The tidal 
stripped material develops a stream forward and backward 
from the satellite along the orbit. In Fig. pj|b, a low surface 
brightness tidal extension of the satellite can be seen at interme- 
diate epochs. At the final snapshot (t = 10.5 Gyr), the residual 
core of the satellite is still detectable, it has ^ 15% of the ini- 
tial mass (Fig. |l2|b) and a central projected surface brightness 
fio ^ 23.8 mag arcsec^^ assuming / L)-^^^^ — 2 M0/L0. 
The initial central surface brightness for the same (M/L)j.g^j 
ratio was 21.3 mag arcsec^^. The stream formed close to the 
satellite has a surface brightness which is 5.0-6.5 mag fainter 
than the center of the satellite. This stream is similar to those 
formed by extra-tidal stars observed close to some DSph satel- 
lites. 

At some snapshots (i.e. 3.7 Gyr, 4.83 Gyr, 6.2 Gyr, 7.3 Gyr, 
8.45 Gyr and 9.6 Gyr), when the tidal streams of the satellite 
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Fig. 9. s-Bl model: a surface brightness fiy assuming 
(7\f/L)j.g^j y = 2 M0/L0, b velocity dispersion of the satel- 
lite as a function of the radius, c distance to the satellite the 
Galactic Center, d evolution of the radii which enclose 70, 60, 
10% of the satellite mass (lines from left to right), and e 
evolution of the mass-to-luminosity ratio measured from the 
velocity dispersion. The observational values of fi and a are 
given in Table |l] 



are akeady formed, the radial velocity, Vr, and the gradient of 
the radial velocity, \dv/dh\, of the s-B2b satellite and of the 
stream around the satellite reproduce the observations of Sgr 
DSph region. The observational values of and \dv/dh\ of 
Sgr are given in Table |l]. However, the values of both quanti- 
ties in the models are strongly dependent on the position and 
orientation of the orbit. Therefore, small perturbations in the 
orbit of the models could lead to differences in Vj. of several 10 
km s^^. In general, the sign of the variation of the radial veloc- 
ity depends on the position of the satellite along the orbit. For 
some particular positions, the radial velocity gradient is almost 
km s^ ^/degree. The s-B2b satellite at time 8.47 Gyr is an ex- 
ample of a model which has a good agreement in position with 
the observations of Sgr DSph (Fig. [TJc). From the kinematical 
point of view, the particles in the stream show a velocity gra- 
dient \dv/dh\ ~ 4 km s^^/degree (Fig. p4^), which is slightly 
higher than the observed value {\dv/dh\ = 3 km s~ ^/degree), 
whereas the radial velocity is slightly lower Moreover, this 
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Fig. 10. Evolution of the log of the ratio between the mass- 
to-luminosity relation, {M/L)o, inferred from the velocity 
dispersion values, and the real mass-to-luminosity relation, 
(Af/L)j.g^j, for S-B2 satellites. The parameter c = \og{rt/ro) 
is the concentration of the model 



satellite model presents a velocity dispersion (cr 6 km/s. Fig. 
p^) lower than the observed value (a ^ 11 A km/s). However, 
the projected velocity dispersion of the stream remains almost 
constant and equal to the velocity dispersion in the inner region 
of the system. Observational data show that regions close to the 
Sagittarius globular clusters (Terzan 8, Terzan 7, Arp 2 ) have a 
similar projected velocity dispersion (Ibata et al. 1997 1. 

The density contour map of the s-B2b model at 8.47 Gyr 
(Fig. [l4|c) resembles the observational map for the region 
aroun d the c enter of the satellite (see for example Fig. 1 in Ibata 
et al. 1997), in spite of the low resolution (few mass points) 
of the simulation. However, the center of our satellite has a 
steeper density profile, the half -brightness radius of the model 
is Tfib — 0.12 kpc and that of Sagittarius DSph is 0.55 kpc. 
Even at the end of the simulation (t — 10.5 Gyr) the satellite 
remains too concentrated (rhi, = 0.17 kpc). This discrepancy 
could mean that either the real Sagittarius DSph has undergone 
a stronger tidal field, which has caused an effective destruc- 
tion of the satellite core, or it has suffered tidal disruption for a 
longer time, or the initial density profile of the Sagittarius satel- 
lite was more extended than that of our model. In order to test 
the latter hypothesis, we have built the s-B2c model, which is 
less concentrated than s-B2b but more massive. s-B2c follows 
the same orbit as s-B2b. 

s-B2c model: 

The s-B2c satellite survives 5.5 Gyr before disruption. The 
effects at perigalacticon (Fig. p^) are stronger than those of 
the other s-B2 models, because of the lower concentration. Low 
surface brightness trailing and leading streams are formed. At 
5.88 Gyr, the residual satellite looks like a long tidal extension 
with a maximum projected surface brightness ji ^ 27.0 mag 
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Fig. 11. Velocity dispersion, a{r), at four snapshots, for s-B2 
satellites 



The line-of-sight velocity dispersion evolves as shown in 
Fig. [rT| c. In the outer region of the satellite it increases up to an 
almost constant velocity dispersion along the orbit. At the last 
snapshot of Fig. |ll]c (5.88 Gyr) the satellite is akeady disrupted 
and the velocity dispersion of the remnant stream has increased 
with respect to the value before destruction. 

The central {M/L)o ratio obtained from Eq. (|l]) depends 
on the position of the satellite along the orbit (Fig. [lO|). The 
most important variations in the {M/L)o ratio are observed at 
the perigalactic passage. At this time the satellite is compressed 
and stretched, its internal distribution is modified and the tidal 
extensions are more pronounced. At the final steps of the evo- 
lution, the more stripped the satellite is, the higher measured 
the {M/L)o is. 

We have analyzed the satellite characteristics at 5.38 Gyr 
from the beginning of the simulation. It corresponds to the last 
approach to the galactic disc and after the last perigalacticon, 
which produces total disruption. At this time, the position (6 = 
-10.7°) and the Galactic velocity (([/, V, W) - (214, :, 170) 
km s~^) of the satellite agree with the observational data. In 
Fig. |l5|, the kinematical distribution (radial velocity and line- 
of-sight velocity dispersion) along the orbit and the contour 
density maps are shown. The half-brightness radius, r^b, and 
central surface brightness, ^o, of the model roughly fit the ob- 
servations (Table |l|). The values we obtain in the s-B2c sim- 
ulation at 5.38 Gyr are rhb ~ 0.4 kpc and ^ 24.3 mag 
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Fig. 12. Radii which enclose 70%, 60%, 50%, 40%, 30%, 20% 
and 10% of the initial mass of the satellite, for s-B2 models 
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Fig. 14. s-B2b satellite, at 8.47 Gyr: a Radial velocity, Vr, along 
the orbit (constant galactic latitude), b velocity dispersion, cr, 
along the orbit and c contour density map at the same scale as 
the a) and b) plots (the star represents the mass center of the 
satellite, located at 6 = —11.1°), the contours correspond to 
^ = 23.2, 23.6, 24.0, 24.7, 26.8 and 28.7 mag arcsec^ from the 
center 




Fig. 13. Surface brightness /iy of the s-B2 satellites, assuming 
W^)real.v = 2Mo/Lo 



arcsec"^, which are slightly lower than observational values, 
but they evolve to larger half -brightness radius and lower sur- 
face brightness as the satellite approaches the Galactic disc 
and gets more disrupted. For example, at t = 5.565 Gyr we 
obtain the best accord between the simulations and the ob- 



servations for rhb and /io (the values of the model are 0.5 
kpc and 25.4 mag arcsec"^, respectively), but the position and 
the velocity do not reproduce the observations (b ^ 20° and 
V, W) ~ (103, 0, 86) km/s). 

The radial velocity gradient along the orbit of s-B2c satel- 
lite depends on the orientation of the trajectory, as in the s-B2b 
model. At 5.38 Gyr we obtain \dv/db\ ^ 1.5 km s~^/degree 
for the trailing stream, which is lower than the observed value. 
However, the radial velocity Vr is in agreement with the obser- 
vations (vr = 171 km s^^). The satellite velocity dispersion 
measured inside 0.5 kpc has decreased (a- 7 km/s), but it is 
still large enough to give [M / L)o ^ 10(M/i)j.g^j. 

It is interesting to analyze not only the region around the 
center of mass of the satellite but also the tidal tails which are 
formed along the orbit. We have done it for several tidal tails 
close to the satellite. In Fig. |l^, we plot the surface brightness 
of the trailing tail of the s-B2c model at various timesteps, as- 
suming (Af/L)j.g^j — 2 M0/LO (a larger value of (-^/i)real 
moves the lines down in the figure and vice versa). The sim- 
ulations are compared with the trailing tail observed by Ma- 



teo et al. (1998) recently. As we can see, except for the region 



closest to the satellite center (< 12°), the simulations repro- 
duce the observations, giving a better result when the satel- 
lite is more disrupted. This supports the hypothesis that Sgr 
DSph is close to its destruction. In our simulations, we also 
obtain a leading tail similar to the trailing one. Unfortunately, 
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Fig. 15. s-B2c satellite, at 5.38 Gyr: a Radial velocity, Vr, along 
the orbit (constant galactic latitude), b velocity dispersion, cr, 
along the orbit and c contour density map (the star represents 
the mass center of the satellite, located at 6 = —10.7°), the 
solid contours correspond to pi — 24.3, 25.0, 25.526.7, 27.5 
mag arcsec" 

sponds to /i = 30 mag arcsec^^ 



the proximity of the Galactic disc to Sgr center makes diffi- 
cult to test this result of symmetric tails due to still poor ob- 
servational data. To give another example of the complexity of 
the satellite structure, we plot in Fig. the contour density 
levels and the kinematic behaviour of a tidal tail, at 5.25 Gyrs 
from the beginning of the simulation. It looks like the Sgr DSph 
from the kinematical point of view, having radial velocity (Fig. 
[l7|a) and velocity dispersion (Fig. [l7|b) similar to Sagittarius 
observed values (vr — 171 km s^^ and a = 11.4 km s^^). The 
shape of the tail (Fig. [l7|c) also looks like a tidally disrupted 
satellite, however, the central surface brightness is /io ^ 29.5 
mag arcsec"^ (fainter than a typical DSph) and the computed 
{M/L)o ratio is 200(Af/i)j.g^j (higher than the estimation for 
DSph galaxies). Nevertheless, this example illustrates, on one 
hand, the difficulty of distinguishing between a tidal tails and a 
DSph galaxy when the latter is close to the disruption and, on 
the other hand, the possibility of identifying objets (globular 
clusters, etc) which have been tidally disrupted from a satellite 
galaxy by measuring their kinematic characteristics. 
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Fig. 16. Surface brightness /iy of the tidal tail close to the Sgr 
satellite. Lines represent the s-B2c model at various timesteps. 
Filled squares are the data (with the errorbars) by Mateo et al. 
(|T99t) 
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Fig. 17. Tidal tail of the s-B2c satellite: a Radial velocity, Vr, 
along the orbit (constant galactic latitude), b velocity disper- 
sion, (T, along the orbit and c contour density map (the star 
represents the mass center of the satellite, located at & = 5.9°) 
and the contours correspond to ^ = 29.6, 29.9, 30.3, 30.8, 32.0 
mag arcsec^ 



6. Discussion 

Summing up the results obtained in our simulations and com- 
paring them with other publi shed w orks, we agree with other 
authors (Velazquez & White [1995^ Ibata & Lewis |1998[ ) that 
Sgr DSph has a short period orbit (T < 1 Gyr). This constraint 



comes from measurement of the position and the radial velocity 
of Sgr DSph in a realistic model of the Galaxy. However, we 
disagree about the DM content. These previous works claim 
that Sgr DSph is a DM dominated satellite, because, contrary 
to us, they could not obtain a low mass satellite which survives 
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by orbiting in the potential of the Galaxy. Therefore, they can- 
not explain the complex evolutionary pattern inferred for the 
Sagittarius DSph, which has suffered a chemical enrichment 
and evolution and the age interval of its globular clusters which 
suggests a long life orbiting in the Galaxy. 



(i) either the satellite has been orbiting for a long time in a 
more external region (where tidal forces are smaller), it suffers 
dynamical friction or/and a deflection with a dense structure of 
the MW and it reaches the present orbit where, eventually, it is 
disrupted, or (ii) it has formed in a tidal tail which contained an 



Ibata & Lewis (1998) have tested low massive models in or- old steflai- population (Kroupa 1998c) 



der to reproduce the Sgr DSph characteristics, reaching an un- 
successful result on the matter and concluding that Sgr DSph 
must have a large DM content. There are two differences be- 
tween their study and ours. Firstly, they use a rigid potential to 
model the Galaxy, so no energy interchange is allowed between 
the satellite and the Galaxy. This restriction could eventually 
prevent a readjusting of the internal energy of the satellite in 
order to reach an equilibrium with the environment. Even if the 
energy transfer from the satellite to the halo is not important in 
all cases, we emphasize that our treatment is a priori adaptable 
to satellit es of various masses. Secondly (and mainly), Ibata 
& Lewis (1998) do not take into account the tidal potential of 
the Galaxy as they build their satellite model. We have con- 
firmed that the fate of such satellites is different from that of 
the satellites in equilibrium with the tidal potential as consid- 
ered here. Consequently, a more accurate model of the initial 
satellite, that reflects the true dynamical situation of the DSph, 
is required in order to avoid spurious effects in the simulations. 
These differences in the models are responsible for the oppo- 
site conclusions reached by Ibata & Lewis, compared with ours 
concerning DM. Furthermore, the recent observations of Sgr 



by Mateo et al. (1998) prove the existence of a long trailing 



tail, which supports the hypothesis that Sgr DSph has suffered 
strong tidal forces, it is close to disruption and, therefore, it is 
not in virial equilibrium. The existence of this tail also reduces 
the inferred {M/L)o ratio. 

Our low DM satellites {modified King's models) establish 
important restrictions on the satellite formation theory. Thus, 
the satellite galaxies grow either in a quasi-isolated region 
(falling slowly on the center of the primary galaxy and having 
time for readjusting their internal structure to the tidal forces), 
or inside the main galaxy, in the tidal tails of a major accretion 
event, which automatically implies equilibrium with the tidal 
forces at the satellite formation epoch. Furthermore, the sur- 
vival time of the equilibrium satellites in the potential of the 
Galaxy depend on the initial concentration of the DSph (the 
larger concentration, the longer the life-time). We have obtain 
models (i.e., s-B2a) which survive more than 10 Gyr The evo- 
lution of these satellites gives rise to tidal streams and modi- 
fications on the outer velocity dispersion of the satellite and it 
leads to high observed mass-to-light ratios if dynamical equi- 
librium is assumed. We have obtained a rather good qualita- 
tive agreement with the observational constraints for the model 
s-B2c in a time interval (from 5.38 to 5.7 Gyr), although this 
agreement is not achieved simultaneously on all the constraints. 
We remind that the simulation time does not necessarily corre- 
spond with the age of the oldest globular cluster of Sgr and, 
therefore, the s-B2c model is not necessarily in disagreement 
with the observations. In this case, two possibilities could be 
invoked to explain the age of the oldest Sgr globular clusters: 



The tidal tail models illustrate that a DSph without DM or- 
biting through central regions of the Galactic potential could 
preserve its structural parameters for a long time. Later on, the 
DSph could show high inferred {M/L)o values as it becomes 
disrupted. The only restriction of the model is that the satel- 
lite must reach the central regions of the Galaxy in equilibrium 
with the tidal forces. That could be possible whether the satel- 
lite has been formed in a tidal tail of an accretion event or it 
has been slowly accreted from the external parts of the Galaxy. 
The latter possibility could increased the estimated life-time of 
the satellite before destruction. 

Recently, K98 has discussed the parameter space of the 
DSphs which could survive in a tidal field. In spite of the dif- 
ferent initial distribution functions, our s-models seem to be 
located in the available region of the (Mgat, Tsat) space for sur- 
viving satellites. However, the main difference between the ap- 
proaches to the problem is that we build the model of the satel- 
lite initially in equilibrium with the tidal forces by solving the 
Boltzmann equation (Gomez-Flechoso & Domfnguez-Tenreiro 



1998) and that assures a priori the survival of the satellite for a 



long time. 
7. Summary 

We have performed simulations of Sagittarius-like satellites to 
follow the evolution of these systems. The satellite and the pri- 
mary galaxy are modelled as N-body systems, that prevents 
possible purely numerical effects which could appear in a non- 
selfgraviting scheme (i.e. if one of the interacting system is rep- 
resented by a rigid potential). Two satellite models have been 
tested: 

1. First of all, we have used isolated King's models (three 
free parameter models) for the satellite, but they are tidally 
stripped and destroyed in a short interval of time, undergo- 
ing a fast evolution. 

2. In order to correctly model the observed situation we had to 
introduce the modified King's models (two free parameters- 
models), which allow us to build up the distribution func- 
tion of the satellite, taking into account the tidal force of 
the environment. These s-models could correspond to no 
DM dominated satellites. In this framework, the initial con- 
centration of the satellite and the trajectory determine its 
life-time. Central circular orbits lead to more effective tidal 
destruction events than highly eccentric orbits. For a given 
kind of orbit, the most concentrated satellites live longer. 
The model which better reproduces the observations is the 
s-B2c model after orbiting 5.38-5.6 Gyr in the potential 
of the Milky Way. For the same orbits, the s-models sur- 
vive a longer time orbiting in the Galactic potential than 
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the f-models, despite the s-models having smaller masses 
(smaller binding energies). 

According to the results presented here, our main conclu- 
sions about the evolution of Sgr are: 

1. Sagittarius may not be a dominated DM satellite. 

2. It follows an eccentric orbit (perigalacticon and apogalacti- 
con are approximately 15 and 70 kpc). 

3. It was more concentrated in the past than at the present 
epoch. 

4. It could have been orbiting in the Galactic halo for a long 
time (minimum 5 Gyr). 

As general conclusions about the DSph satellites, we can 
remark that: 

1. High values of {M/L)o found in the literat ure (Ib ata et al. 

Mateo |1998| for a re- 



1995; Irwin & Hatzidimitriou 1995 



view of the Galactic satellites) could be due to an erroneous 
use of the virial theorem and not to the presence of DM (as 
already suggested by K97; KK98). 

2. The progenitor of these satellite could be either lumps 
formed in a major accretion event or other dwarf galax- 
ies (maybe dwarf irregulars or dwarf spirals) which fall in 
the potential well of the Galaxy, lose their gas content and 
evolve to the equilibrium configuration in the tidal field of 
the MiUcy Way, giving rise to the population of DSphs. The 
finding of intermediate states of this evolution could be a 
support to the proposed scenario. 

3. According to some of the present results, some DSphs are 
able to survive for long times orbiting in the Galactic halo. 
That should be taken into account for the calculations of the 
dissolution rate of dwarf galaxies in the halo of the primary 
galaxy. 
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